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ABSTRACT 

We generate synthetic Hi Galactic plane surveys from spiral galaxy simulations 
which include stellar feedback processes. Compared to a model without feedback we 
find an increased scale height of Hi emission (in better agreement with observations) 
i-£h and more realistic spatial structure (including supernova blown bubbles). The syn- 

thetic data show Hi self-absorption with a morphology similar to that seen in observa- 
Q tions. The density and temperature of the material responsible for Hi self-absorption 

?-) is consistent with observationally determined values, and is found to be only weakly 

c/3 dependent on absorption strength and star formation efficiency. 
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1 INTRODUCTION 

Hi emission, from the 21-cm line of atomic hydrogen, is 
widely used as a tracer of large scale Galactic structure 
( Kalberla & Kerp 2009 ) but also reveals smaller scale struc- 
ture in the interstellar medium. This smaller scale struc- 
ture includes shells formed by stellar feedback processes, 
which inject kinetic and thermal energy, along with heavy 
elements. Such feedback is instrumental in determining the 
evolution of a galaxy and its effects can be seen in our own 



Galaxy (Heiles 1979 1984 Taylor et al. 20031 and also in 



nearby external galaxies ( Fukui et al.|2 009 Bagctakos et al 



|2011[ ). Hence feedback processes are seen to play an im- 
portant role in the formation of smaller scale structure and 
its appearance in Hi emission. Numerical simulations of the 
interstellar medium (ISM) in disk galaxies are becoming in- 
creasingly sophisticated and can now include the effects of 



stellar feedback processes (Wada 2008 Shetty & Ostriker 



[2008{|Dobbs et al.|20lTb||Tasker fc Tan|2011[ ). At the same 
time advances are being made in the generation of syn- 



thetic observations from simulations ( |Hennebelle et al. 2007 
|Douglas et al.|[2010| |Parkin|[2011| and references therein). 
These synthetic observables allow a direct comparison be- 
tween simulation results and real observations, and as the 
models become increasingly sophisticated we can expect the 
synthetic observations to compare more favourably with real 
observations. As well as testing whether models can produce 
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realistic observables, such synthetic observations also allow 
features to be related more directly to underlying physical 
processes than is possible with real observations. Hence syn- 
thetic observations can be used to test physical mechanisms 
proposed to explain features seen in observations. 

The transition of hydrogen gas from the atomic phase 
to the molecular phase is a key step to molecular cloud, and 
ultimately star formation. In order for this transition to oc- 
cur the gas must become cooler and denser than typical ISM 
material. The passage of material through the shocks associ- 
ated with the density waves in a grand design spiral galaxy 
results in substantial compression ( |Roberts|1972| , which al- 
lows the gas to cool efficiently down to low temperatures 



(|Cowie|1981||Bergin et al.|2004||Dobbs et al.|2008"}|Kim et~al 
|2008[ ). The high density, low temperature material in a spi- 
ral shock provides favourable conditions for the formation of 
molecular clouds (although molecular clouds can also form 
in the absence of spiral shocks). Whilst CO traces molec- 
ular clouds, cold Hi with a low molecular fraction is more 
difficult to detect. It can however be detected as Hi self- 
absorption (HIS A). HIS A is observed when colder, denser 
foreground material is located along the same line of sight, 
and at the same line of sight velocity, as warmer background 
material ( |Li fc Goldsmith||2003| |Gibson et al]|2005[ ). Given 
these requirements, using HISA is not a suitable method to 
obtain a complete map of cold Hi. However HISA can pro- 
vide information in cases where other tracers of star forming 
clouds (such as the J=l-0 rotational transition of the CO 
molecule) are not effective ( Douglas fc Taylor|2007 l, for ex- 
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ample shortly after passing through a shock (Bergin et al 



2004 1, or in photon dominated regions where there is less ef- 



fective self-shielding of CO against the interstellar radiation 
field compared to H2 ( jKaufman et al.|1999|| Allen et al. 2004 
Shetty et aL]|2011[ ). Thus although HIS A i s found to corre- 
late with emission from molecular species ( Goldsmith & Li 
2005 Kavars et al.|200"5 1 the correspondence between HIS A 
and tracers of star forming clouds is not straightforward. 
Synthetic observations provide a valuable way of examining 
tracers such as HIS A in a more controlled environment. 

In this paper we present synthetic Hi Galactic plane 
surveys derived from simulations of galaxies which include 
stellar feedback and compare the Hi emission and HISA 
structure with similar features found in observations. We 
begin in Section [2] by outlining the method used to make 
the synthetic Hi surveys. In Section[3]we assess the effect of 
feedback on the observable properties of the Hi gas by com- 
paring our synthetic data with observational data and with 
a previous synthetic survey which did not include feedback. 
The distribution of HISA in our synthetic observations is 
discussed in Section |4.1| and the properties of the material 
producing HISA are investigated in Section |4.2| We finish 
by presenting our conclusions in Section [5] 



2 METHOD 

Our synthetic observations are generated from an SPH 
(smoothed particle hydrodynamics) simulation of a spiral 
galaxy. The variable spatial resolution of the SPH method 
permits a large range of density and spatial scales to be re- 
solved. This enables a simulation with a domain covering 
a whole galaxy, which can still resolve individual molecular 
clouds. SPH is a Lagrangian method which means that the 
material associated with an SPH particle retains its identity 
as the simulation progresses. Consequently we can track ma- 
terial as it passes thorough a spiral arm and determine its 
fate at a later time. A synthetic Galactic plane survey was 
previously presented by |Douglas et al!] (2010) (hereafter Pa- 
per 1) using an SPH model without feedback processes. The 
new results presented here include stellar feedback processes 
in the SPH model, which was proposed as an important 
mechanism to resolve differences between Paper 1 results 
and observations. 

The SPH simulations used for this work are models C 



and D from Dobbs et al. (2011b I which include stellar feed- 



back and self- gravity of the gas, which were not included in 



the simulation of Dobbs (20081 used in Paper 1. In all the 



ergy, where gas is assumed to have formed stars, although 
sink or star particles are not introduced. For feedback to oc- 
cur the density of a gas particle in a converging flow must ex- 
ceed 1000 cm -3 , whilst the surrounding region of gas (within 
a radius of ~ 20 pc) must be gravitationally bound and in 
an energetically favourable state ( Dobbs et al.pOlla l. The 
number of supernovae is estimated from the mass of gas 
in the computed region, adopting a certain star formation 
efficiency and a Salpeter initial mass function (IMF). For 
each feedback event, the total amount of energy deposited 



is given by equation 1 from Dobbs et al. (2011al, which 



computes the number of massive stars expected to form, 
and assumes each supernova injects 10 51 erg of energy. The 
number of stars formed (and therefore the amount of energy 
deposited) depends on a star formation efficiency parameter, 
e, which represents the fraction of the molecular mass in each 
bound region assumed to be turned into stars. The star for- 
mation efficiency is an absolute value and does not depend 
on the free-fall time. The energy is distributed according to 
a snowplough solution, which describes the pressure-driven 



phase of a supernova after the blast wave (Woltjer 1972 



Chevalier 1974 McKee & Ostriker 19771 and is deposited 



as 2/3 kinetic energy, and 1/3 thermal energy. The feedback 
is assumed to be instantaneous, and although we consider 
supernova explosions, the energy could account for numer- 
ous feedback processes, such as stellar winds, radiation and 
supernovae. 

Model data at a simulated time of 250 Myr are used, 
at which time the simulated galaxy has reached a state of 
quasi-equilibrium. In Model C approximately one third of 
the ISM is in each of the cold (< 150K), unstable (150- 
5000K) and warm (> 5000K) phases. With a higher star for- 
mation efficiency there is more material in the warm phase 
and less material in the cold phase. The thermal evolution 
of the gas is modelled according to the thermodynamics of 



Dobbs et al. (20081. The gas exhibits a maximum temper- 



ature of 2 x 10°K. A lower limit of 20K is imposed prior 
to the temperature update in the model. This limit pre- 
vents material from becoming too cold for the simulation to 
treat accurately but allows some cooling below 20K to oc- 
cur. The simulation follows the evolution of molecular gas, 
although only a small fraction (< 10 percent) is molecular. 
Energy is injected as soon as the criteria for feedback are 
met, which results in dense gas being rapidly disrupted. In- 
cluding a delay before injecting energy from feedback would 
enable longer for molecular gas to form. 

The SPH data provide three dimensional distributions 
of Hi density, temperature and velocity which are used to set 
up an adaptive mesh refinement (AMR) grid for the TORUS 



radiative transfer code (Harries 2000). TORUS carries out 



simulations the gas is assumed to orbit in a fixed potential, 
composed of halo, disc and four armed spiral components. 
The resulting rotation curve is comparable to that of the 
Milky Way, which allows structures in the simulations to be 
compared to observed structures in our own Galaxy. For ex- 
ample the simulations have spiral arms which are similar to 
those in our Galaxy which allows us to generate synthetic 
images analogous to observations of the nearby Perseus arm 
in our Galaxy. The total surface density (including Helium) 
is constant with values of 8 Mqpc -2 for models including 

feedback and 10 Mqpc -2 for the model without feedback et al. (20101 using the method of Rundle et al. (20101. The 
used in Paper 1. This is comparable to the surface density 
in the solar neig hbourhood ( jWolfire et al.|2003[ ). 

Stellar feedback is inserted as kinetic and thermal en- 



non-parallel ray traces (one per image pixel and velocity 
channel) to solve the radiative transfer equation and gener- 
ate three-dimensional spectral data cubes of 21 cm Hi emis- 
sion in Galactic longitude-latitude-velocity co-ordinates. 

The AMR grid is generated by splitting the grid cells 
so that there is never more than one SPH particle per cell. 
Mapping temperature, density and velocity onto the AMR 
grid is carried out as described in Paper 1 an d |Acreman| 



observer is placed within the model galaxy, at a location 
analogous to that of the Sun within our own Galaxy, and 
the synthetic survey is generated as if it were the Galactic 
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second quadrant (90° < I < 180°). The ray tracing uses 



the density sub-sampling method of Rundle et al. (20101 to 



linearly interpolate density values within a cell on the AMR 
grid. Density sub-sampling was not required in Paper 1 as 
the spatial resolution of the SPH simulation and AMR grid 
were finer (the computational demands of adding feedback 
processes require that a lower spatial resolution is used). 
The data cubes have a velocity resolution of 0.5 km/s and 
an angular resolution of 1 arcmin in longitude and latitude. 

Each cell on the AMR grid is assigned a thermal line 
width, which does not include any additional turbulent com- 
ponent to account for unresolved structure. This avoids 
adding an ad-hoc parameter but will tend to under-estimate 
the line width if unresolved structure makes a significant 



contribution to the velocity dispersion (Hennebelle et al 



2007 1. The velocity dispersion of the gas and clouds in the 



simulations are discussed in Dobbs et al. (2011b) 



3 EFFECT OF FEEDBACK ON HI 
STRUCTURE 

3.1 Latitude-longitude structure 

Figure [T] shows latitude-longitude plots of Hi emission for 
simulations with and without feedback, and also for the 



Canadian Galactic Plane Survey (CGPS) observations (Tay 
lor et alT||2003 > . The CGPS data are from a velocity chan- 
nel which contains emission from Perseus arm material. In 
the synthetic observations the observer has been positioned 
so that at similar velocities we see a spiral arm analogous 
to the Perseus arm. Fig. |l(a)| shows synthetic observations 
from a simulation without feedback from Paper 1 (hereafter 



referred to as NoFeedback). Figure 1(b) is from Run C of 
|Dobbs et al.| ( [20"llb| and has feedback included with 5% star 
formation efficiency (hereafter referred to as Feedback5), 
and Fig. 1(c) is from Run D of Dobbs et al. (2011b I and 



has feedback included with 10% star formation efficiency 
(hereafter referred to as FeedbacklO). The star formation 
efficiency referred to here is the e parameter described in 
Section |1 

The CGPS data show a broad region of emission which, 
in places, extends beyond ±2 degrees outside the mid-plane. 
In contrast the NoFeedback run has emission which is much 
more confined to the mid-plane, with bright ridges of emis- 
sion at ±1 degree which are not seen in the CGPS data. 
Without feedback the gas in the model galaxy is overly con- 
fined to the mid-plane, resulting in large mid-plane opti- 
cal depths. In the NoFeedback model the column density 
is highest within ±1 degree of the mid-plane but the accu- 
mulation of cold, dense material results in excessively high 
absorption and correspondingly low levels of Hi emission (in 
the optically thick limit the brightness temperature will sat- 
urate at the spin temperature). When feedback is included 
(Fig. 1(b) and 1(c) I the gas in the model galaxy is much less 
confined to the mid-plane and the bright ridges of emission, 



seen in Fig. 1(a) at about 1 degree above and below the mid- 



plane, are not present. The Hi emission now extends further 
out of the mid-plane, in better agreement with the CGPS 
observations. 

Feedback results in significant holes in Hi emission, with 
FeedbacklO in particular having much less contiguous emis- 
sion than the CGPS observations. The CGPS observations 



have small scale filamentary structure, which is not seen in 
the synthetic observations, however, we do not expect to re- 
produce this structure at the current SPH resolution. The 
resolution of the SPH simulation is governed by the smooth- 
ing lengths of the particle^] The density threshold for feed- 
back to occur is at a number density of 1000cm -3 which 
corresponds to a smoothing length of 5.6 pc (with a particle 
mass of 25OOM0). At a distance of 2.5 kpc (typical of mate- 
rial seen in Fig[l]) this corresponds to an angular size of 0.13 
degrees (or 7.8 arcmin). Structure on this scale is seen in the 
synthetic data (e.g. at I = 120 in Fig 1(c) I but is spherical, 
rather than filamentary. 

3.2 Vertical distribution of Hi emission 

In order to allow a quantitative comparison of the vertical 
distribution of Hi emission, longitudinally averaged profiles 
of brightness temperature against latitude were extracted for 
the longitude range 126-144 degrees, in the same velocity 
channels shown in Fig. [I] Emission seen in this longitude 
and velocity range is from spiral arm material at a distance 
of approximately 2.5 kpc from the observer. At this distance 
material at a latitude of 1 degree is 44 pc out of the mid- 
plane. 

The spiral arms in the simulations have a narrow ve- 
locity width, compared to the observed Perseus arm, and 
this results in emission from the simulated spiral arms 
being distributed over fewer velocity channels but with 
increased brightness temperatures (the structure of the 
arms in longitude-velocity space is discussed further in Sec- 
tion 3.31. Consequently the synthetic surveys have higher 
brightness temperatures than the CGPS data and the shape 
of the raw profiles cannot easily be compared. To facilitate 
a quantitative comparison the synthetic profiles were scaled 
by a constant scaling factor chosen to minimise the RMS 
difference from the CGPS data. The expression used to find 
the best value of the scaling factor / was 



RMS = 



1 N 



N 



(1) 



where TV is the number of points in the CGPS profile, d 
is the brightness temperature in bin i of the CGPS profile, 
and Sj is the brightness temperature in bin j of the synthetic 
profile, where bin j and bin i correspond in latitude. Every 
point in the brightness temperature profile was then scaled 
by the value of / which minimised the expression in Eqn. [I] 
As the profiles are asymmetric the fit was repeated with 
the latitude axis inverted. The scaling factors and resultant 
RMS differences from the CGPS profile are shown in Table[l] 
The scaled profiles and the unmodified CGPS profile are 
plotted in Fig. |2(a)| and in Fig. |2(b)"] (with a reversed latitude 
axis in the second plot). 

The fit for both the feedback runs is significantly bet- 
ter when the latitude axis is reversed, as the skewness of 

1 The smoothing length h is given by 

. V3 



where m is the particle mass, p is the particle density, and r\ = 1.2 
( Price fc Monaghan|2007 l. 
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Figure 1. Brightness temperature of Hi emission, in Galactic latitude and longitude co-ordinates, from the Perseus arm in the CGPS and 
Perseus arm analogues in the synthetic data. The longitude coverage is the same in each case but the latitude coverage varies according 
to the extent of Hi emission (synthetic data) or survey coverage (CGPS). Fig. |l(a)]is fr om a model galaxy without feedback, Fig. |l(b)] is 
from a model galaxy with feedback and 5 per cent star formation efficiency, Fig. |l(c)| is from a model galaxy with feedback and 10 per 
cent star formation efficiency, and Fig. 1 1(d)] is from the Canadian Galactic Plane Survey. 
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(a) Latitude axis not inverted 



(b) Latitude axis inverted 



Figure 2. Profiles of Hi brightness temperature with latitude. The solid line is from the CGPS, the dashed line is from the NoFeedback 
model, the dotted line is from Feedback5 and the dot-dashed line is from FeedbacklO. The normalisation of the synthetic profiles has 
been scaled to fit the CGPS data, so that the shape of the profiles can be readily compared. In Fig. |2(b)| the latitude axis is inverted so 
that the synthetic profiles better match skewness of the CGPS profiles. 



Table 1. Scaling factors applied to latitude profiles from syn- 
thetic data and the corresponding RMS differences from the 
CGPS profile. 

Simulation Scale factor RMS Reversed axis 

difference (K) axis? 



NoFeedback 


0.82 


28.6 


No 


Feedback5 


0.535 


8.86 


No 


FeedbacklO 


0.762 


8.20 


No 


NoFeedback 


0.821 


28.5 


Yes 


Feedback5 


0.539 


3.39 


Yes 


FeedbacklO 


0.767 


3.51 


Yes 



the profile better matches the skewness of the CGPS data. 
Both Feedback^ and FeedbacklO fit the CGPS data well, 
provided the normalisation of the profile is scaled, and it is 
not possible to infer from these results whether 5% or 10% 
efficiency is a better fit to the data. Conversely the NoFeed- 
back model does not fit the CGPS data well as the profile 
shows a large dip around the mid-plane giving a qualitatively 
different profile to that observed. 

Although the star formation efficiency is different in 
Feedback5 and FeedbacklO, the rate of star formation is 
comparable at 250 Myr, as the process is largely self- 
regulating (as shown in fig. 11 of |Dobbs et al.|2011a[ ). Conse- 
quently the energy input from stellar feedback is similar for 
Feedback5 and FeedbacklO (as the supernova rate is propor- 
tional to the star formation rate) and we expect the impact 
on the vertical distribution of Hi to be similar. Conversely 
a synthetic survey from a simulation with twice the surface 
density (run M10 from |Dobbs et al.|2011a| not shown here) 
has a larger scale height than Feedback5 and FeedbacklO. 
This is consistent with the higher star formation rate (ap- 
proximately four times higher) and an increased level of stel- 
lar feedback. 



The amount of energy injected into the ISM in the sim- 
ulations is directly proportional to the star formation rate. 
The star formation rate (SFR) in the Feedback5 and Feed- 
backlO models is ~ O.7M0 /year, which is consistent with the 
lower limit of the observationally determined Galactic star 



formation rate of Robitaille & Whitney 



(2010) but approx- 



Murray & Ra hman] 



imately a factor 2 less than the SFR of 
(20101. We note that SFRs determined from both obser- 
vations and models are sensitive to assumptions about the 
slope and limits of the IMF ( jCalzetti et al.||2009[ ). Addi- 
tionally we can only insert feedback on a resolvable scale in 
the simulations, and cannot take into account all the sub- 
resolution processes which may dissipate or radiate energy. 
Hence we cannot definitively say how well energy injection 
in the models is similar to that in the Milky Way. However 
the effect of feedback has clearly been to make the models 
more realistic in this regard. 



3.3 Longitude-velocity structure 

Longitude-velocity plots in the Galactic mid-plane are 
shown in Fig. [3] These plots show local material around 
v=0 and Perseus arm material at more negative velocities 
(between —20 and —60 km/s). The inclusion of feedback re- 
sults in much more structure in the spiral arms in longitude- 
velocity space. In Feedback5 the Perseus arm material is con- 
tiguous but structured, whereas in FeedbacklO the Perseus 
arm material appears less contiguous. The local material is 
highly disrupted by feedback but as this material is close to 
the observer the angular size of the SPH smoothing length 
is large and the material is poorly resolved. 

The CGPS data show material at more negative veloc- 
ities than the simulations (beyond —60 km/s) but we would 
not expect to see material in this velocity range in the syn- 
thetic data because the model galaxies do not have a suffi- 
ciently large radial extent to represent this material. For gas 
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Figure 3. Brightness temperature of Hi emission, in Galactic longitude- velocity co-ordinates in the Galactic mid-plane. Fig. |3(a)] is from 
a model galaxy without feedback, Fig. |3(b)| is from a model galaxy with feedback and 5 per cent star formation efficiency, Fig. |3(c)| is 
from a model galaxy with feedback and 10 per cent star formation efficiency, and Fig. |3(d)]is from the Canadian Galactic Plane Survey. 
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in axisymmetric circular rotation the line of sight velocity 



(from equation f of Kalberla & Kerp (20091) is 
^0 , 



v(R,z) = 



R 



■9(7?, 



sin I cos b 



(2) 



where v(R,z) is the line of sight velocity at point (R,z), 
O (R, z) is the tangential velocity at (R,z), 0© is the ob- 
server's tangential velocity, Rq is the distance of the ob- 
server from the Galactic centre, and I and b are Galactic 
latitude and longitude. In our synthetic surveys the observer 
is located at Rq = 7.1 kpc within a model galaxy with an 
outer extent of 10 kpc, and at both these locations the tan- 
gential velocity is close to 220 km/s. Hence for material in 
the Galactic plane (b — 0) the line of sight velocity at the 
outer edge of the model galaxy is 



v (R, z) = — 63.8sinZ km/s 



(3) 



At / = 110 the maximum extent in velocity space is expected 
to be —60 km/s decreasing to —22 km/s at I — 160. This is 
in agreement with the most extreme negative velocities seen 
in the synthetic data in Fig. [3] 



4 HI SELF-ABSORPTION 

When generating synthetic observations the calculation of 
Hi intensity can be split into separate emitting and absorb- 
ing components. This allows us to produce a data cube con- 
taining only the absorption component, in order to identify 
where HISA is present. Furthermore, each cell of the AMR 
grid is identified as either a net source of emission or a net 
source of absorption by calculating the change in intensity 
due to the grid cell, normalised by the column density of the 
cell. The SPH particles are then associated with the change 
in emission from the AMR cell in which they reside, thus we 
are able to determine which SPH particles are responsible 
for HISA. Compared to observers we are in the privileged 
position of being able to unambiguously identify absorbing 
components in the data cube and furthermore we can iden- 
tify absorption with specific SPH particles. In Section |4.1| wc 
use the absorption-only data cubes to examine the distribu- 
tion of HISA on the plane of the sky, then in Section [4. 2| we 
use the identification of absorbing SPH particles to investi- 
gate the properties of the material which causes HISA. 



3.4 Expanding shells 

Several shells of material are seen in the synthetic obser- 
vations (see Fig. [4] for an example from the Feedback5 
model), similar to shells seen in Hi observations ( |Heiles|1979| 
1984) |Hu||1981| |McClure-GrifEths et al^[2^02^|FdTlerova fc| 
Palous|2005 1 . Such structures are expected as a consequence 
of energy feedback from massive stars and the feature in 
Fig. [4] is associated with a feedback event which occurred 
1.2 Myr earliei^] The feedback event comprised 20 super- 
novae, an atypically energetic event, which injected a total 
of 2 x 10 52 erg of energy into the ISM. 



4.1 HISA distribution 



In longitude-latitude space (Fig. 4(a) I an expanding 



shell appears as a variable radius ring with a maximum ra- 
dius at the central velocity of the material. The radius of the 
shell decreases in velocity channels away from the central 
value and terminates with two "caps" of emission from the 
front and back of the bubble. By calculating the average sur- 
face brightness in concentric annuli about the centre of the 
shell, the structure can be plotted in velocity-radius space. 
A plot of this type, made using the ksh ell to ol from the 
KARMA software package]^] is shown in Fig. 4(b) In velocity- 
radius space the structure is an arch shape which confirms 
that this structure is indeed expanding. 



Figure 4(b) shows that the shell has an expansion veloc- 



ity of approximately 10 km/s and a maximum angular radius 
of approximately 2 degrees, which corresponds to 80 pc at 
2.3 kpc (the distance from the observer's position). The shell 
radius and expansion velocity are smaller than typical values 



from the samples of Heiles ( 1979) and Heiles ( 1984 1, however 



Heiles ( 1984 1 notes a bias towards selecting larger shells and 



our shell is more typical of the shells found by Hu ( 1981 1 



The size of our shell is also consistent with the smaller shells 



in the sample of McClure-Griffiths et al. (2002), observed in 



the Southern Galactic Plane Survey. 

2 This time does not include the age of the supernova bubble 
when the energy is inserted. This timescale (denoted t in Ap- 



pendix 1, Dobbs et al.|2011a l is of order 10 5 years. 
3 http:/ /www. science-sottware.net/karma/ 



The velocity integrated absorption component is plotted in 
Fig.[5]for the NoFeedback, Feedback5 and FeedbacklO mod- 
els. Equivalent plots from the CGPS data are shown in fig- 
ure 1 of |Gibson et al.| ( |2005 l. 



In the NoFeedback model (Fig. 5(a) ) there is a broad 
band of strong HISA, due to an over concentration of Hi 
in the mid-plane. HISA in the NoFeedback case is confined 
to ±1 degree from the mid-plane, whereas observed HISA 
has a greater vertical extent. The HISA morphology in the 
Feedback5 model is more realistic, with extended, low inten- 
sity HISA surrounding knots of stronger absorption. Both 
the morphology and vertical extent of the Feedback5 model 
are more similar to the CGPS observations of Gibson et al. 
( |2005 1. Although the Feedback5 HISA has less substructure 
in the strong HISA complexes this is to be expected, given 
the small spatial scale of the observed structures relative to 
the model resolution. In the FeedbacklO model the HISA has 
a much reduced diffuse component, compared to the other 



models. Gibson et al. (2005) find nearly ubiquitous weak 



HISA and in this regard the FeedbacklO model does not 
match the observations. Models with higher star formation 
efficiencies have less material in the cold phase, as shown in 



fig. 4 of Dobbs et al. (2011al, so FeedbacklO is expected to 



have less HISA than Feedback5 as there is less cold atomic 
hydrogen. We conclude that the Feedback5 model has a more 
realistic HISA morphology than either of the other models, 
albeit with limited spatial resolution. 



4.2 Properties of material responsible for HISA 



Figure [6] shows a histogram of number densities (Fig. |6(a)| 
and temperatures (Fig. 6(b) I for absorbing particles from 



Feedback5, located in the region used to generate the syn- 
thetic survey. The solid line is for all particles in a cell 
with net absorption (6% of particles in region), the dashed 
line is for particles associated with absorption stronger than 
1 x 10" 23 erg/s/sr (2% of particles in region) and the dot- 
ted line is for particles associated with absorption stronger 
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Galactic longitude (degrees) Radius (degrees) 

(a) Shell in longitude-latitude space (b) Shell in velocity-radius space 

Figure 4. An expanding shell of Hi emission from the Feedback5 model. The velocity-radius plot (right) shows the arch shape which 
is expected from an expanding shell. The angular offset is the radius of an annulus about the centre of the cell, where the centre is 
determined by eye. 



than 2 x 10 -23 erg/s/sr (1% of particles in region). There 
is no evidence for particles associated with stronger HISA 
to be higher density or at a lower temperature, indeed the 
temperature range is 10-100K (typical of the cold neutral 
interstellar medium) for almost all the absorbing particles. 
For these histograms the mean Hi number density is 130-190 
cm -3 and the mean temperature is 32-36 K. 

The density and temperature distributions for absorb- 
ing particles from FeedbacklO are shown in Fig. 6(c) and 
Fig. |6(d)| respectively. Particles associated with net absorp- 
tion are plotted as a solid line but no other thresholds are 
used, as the number of absorbing particles is much smaller 
than for Feedbacks (1% of particles are absorbing in Feed- 
backlO). For FeedbacklO the mean Hi number density is 



there is sufficient spread in our temperature values to en- 
compass their observationally determined temperatures (see 



Fig. 6(b) and Fig. 6(d) I. Material can drop out of the gas 
phase and form supernovae above a density of 1000 cm -3 , so 
the densest and coldest regions of HISA producing gas (seen 
in the observations) are at the limit of the SPH simulations' 
representation. 

The next point to be addressed is whether particles 
which contribute to HISA are in a different phase to other 
particles at similar temperatures, and whether HISA pro- 
ducing material is in pressure balance with other material. 
In Fig. [7] we plot pressure against number density for HISA 
particles (Fig. |7(a')j ) and for all other particles with a temper- 
ature below 150K (Fig. 7(b) I. There are approximately five 

indicat- 



490 cm and the mean temperature is 27K. The range of times as many particles in Fig. 7(b) as in Fig. 7(a 
temperatures and densities seen in FeedbacklO are similar to 
those seen in Feedback5, although the HISA in FeedbacklO 
is from slightly colder and denser material. The star forma- 
tion efficiency has only a modest effect on the properties of 
material seen in HISA, even though the morphology of the 
HISA is significantly different. 

I n thei r study of HISA clouds in Perseus, |Klaassen"| 
et al. ( 2005 I found number densities between 100 cm -3 and 
, and spin temperatures in the range 12K-24K, 



1200 cm" 

in an extended HISA cloud which they term the "complex" . 
In a smaller, isolated HISA feature, which they term the 
"globule" , the spin temperatures were in the range 8K-22K 
(no density determination was made for the globule.) The 
globule is a compact structure (unresolved in the 1 arcmin 
main beam of the CGPS) and is much smaller than the res- 
olution limit of our synthetic observations (see section pTTj ) . 
The complex region is much larger and hence is more like 
the HISA clouds seen in our synthetic data. The density 
of HISA producing material in Feedback5 and FeedbacklO 
is consistent with the density of the complex region found 



by Klaassen et al. (20051. Although the mean temperature 



of our HISA producing material is slightly higher than the 



temperatures found by Klaassen et al. (2005) we note that 



ing that there is a significant amount of cold, dense Hi which 
is not observed in HISA. For particles with a number den- 
sity above 10 2 cm -3 there is no apparent difference between 
the two populations. At lower densities the distributions are 
also similar but there is more spread (towards higher pres- 
sures) in the non-HISA particles. The thermal pressure of 
HISA material does not appear to be significantly different 
to other material, although we note that the effect of velocity 
dispersion has not been included here, which would provide 
support against gravity in addition to thermal pressure. 

Although the material responsible for HISA is identifi- 
able as the cold neutral medium, which is expected to go 
on to form stars, we can be more specific about the fate of 
the SPH particles responsible for HISA. Particles which are 
in giant molecular clouds (GMCs) are identified, using the 
clump finding algorithm of |Dobbs et al.| ( |2008[ ), and corre- 
lated with whether the material is also observed in HISA. 
In the Feedback5 model 4.0% of the HISA particles are 
in GMCs, compared to 2.2% of material in GMCs for the 
galaxy as a whole, and 70% of HISA material is involved 
in a feedback event (i.e. involved in star formation) within 
the next 20 Myr. The fraction of HISA material in clouds, 
and the fraction which forms stars, is not found to vary with 
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Figure 5. Hi self-absorption amplitude, calculated from the synthetic observations, where the absorption is integrated over all velocity 
channels. Fig.|5(a)]is from a model galaxy without feedback, Fig.|5(b)]is from a model galaxy with feedback and 5 per cent star formation 



efficiency, Fig. |5(c)]is from a model galaxy with feedback and 10 per cent star formation efficiency 
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Figure 6. Number density and temperature distributions of SPH particles associated with Hi self-absorption (HIS A). For the Feedback5 
model three different absorption thresholds are used to define HISA particles: dl < (solid line), dl < — f X 10~ 23 crg/s/sr (dashed 
line), and dl < 2 X 10~ 23 erg/s/sr (dotted line). For the FeedbacklO model there are fewer HISA particles and only one line is plotted 
(particles with dl < 0). 



HISA strength. For FccdbacklO we find that 2.5% of HISA 
material is in GMCs, compared to 0.9% in the galaxy as a 
whole, and 57% of HISA material goes on to be involved in 
star formation. 



5 CONCLUSIONS AND DISCUSSION 

The inclusion of feedback in the galaxy models has a sig- 
nificant effect on the derived synthetic Hi Galactic plane 
surveys. In the model without feedback (NoFeeback) gas is 
overly confined to the mid-plane, which results in excessive 
absorption. Consequently there are bands of bright emis- 
sion above and below the mid-plane (which are not seen in 
observations) and the vertical extent of the Hi emission is 
too small. The two models which include feedback (Feed- 
backs and FeedbacklO) both have a larger vertical extent of 



Hi emission and their profiles of Hi emission with latitude 
match the CGPS observations well (provided the normali- 
sation is allowed to vary) . Based on the vertical distribution 
of Hi emission alone it is not possible to determine whether 
Feedback5 or FeedbacklO is more realistic. 

When feedback is included more structure is seen in 
Hi emission, including bubbles of material associated with 
supernova events, however the model does not yet have suf- 
ficient spatial resolution to capture the very fine scale struc- 
ture seen in the CGPS data. With increases in available com- 
puting power it will be feasible to run the SPH simulations 
using more particles, giving a higher spatial resolution. This 
should result in more realistic small scale structure in the 
synthetic observations. However if the filamentary structure 
is due to the presence of magnetic fields then these too will 
need to be included before a good match with the observed 
Hi morphology on a small scale can be expected. 
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Figure 7. Plots of pressure against number density for particles producing Hi self-absorption (Fig. |7(a)| and all other particles with 
temperatures below 150K (Fig, |7(b)| l, showing that HISA particles are not over-pressured compared to other cold particles. 



One further respect in which the synthetic data are 
more realistic is the properties of HISA. The NoFeedback 
model has a broad band of strong HISA, unlike HISA seen 
in observations. When feedback is included with 5 per cent 
star formation efficiency the HISA structure is much more 
like that in the observations of |Gibson et al.| ( |2005[ ) . How- 
ever if the star formation efficiency is increased to 10 per cent 
there is insufficient HISA, in particular there is no ubiqui- 
tous weak HISA as found by |Gibson et al.| ( |2005[ ). Based on 
the HISA morphology we conclude that a 5 per cent star 
formation efficiency is more realistic than a 10 per cent star 
formation efficiency. 

With the simulations with feedback the HISA distri- 
bution is more realistic and it is meaningful to examine 
the properties of material which is seen in HISA. The tem- 
perature and density properties are confined to a well de- 
fined range (typical of cold Hi) which is independent of how 
strongly the material is absorbing, and only weakly depen- 
dent on the star formation efficiency. Although the corre- 
spondence between HISA and molecular clouds is not one- 
to-one it is clear that HISA preferentially selects material in 
GMCs. However the nature of clouds is dynamic and a given 
cloud will not be composed of the same material as it evolves. 
Consequently material which is seen in HISA at any given 
time may become part of a GMC later, and HISA material 
currently in a GMC can leave the cloud at a later time. Un- 
derstanding the complex relationships between HISA, GMC 
material and star formation could be facilitated with a time 
series of synthetic observations, which would track particles, 
determine when they are observable as HISA, and correlate 
this with the molecular gas fraction. 
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